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We provide an algebraic procedure to find the eigenstates of two-charged particles in an oscillator 
potential, known as Hooke's atom. For the planar Hooke's atom, the exact eigenstates and single 
particle densities for arbitrary azimuthal quantum number, are obtained analytically. Information 
entropies associated with the wave functions for the relative motion are then studied systematically, 
since the same incorporates the effect of the Coulomb interaction. The quantum pottery of the 
information entropy density reveals a number of intricate structures, which differ significantly for 
the attractive and repulsive cases. We indicate the procedure to obtain the approximate eigen states. 
Making use of the relationship of this dynamical system with the quasi-exactly solvable systems, 
one can also develop a suitable perturbation theory, involving the Coulomb coupling Z, for the 
' approximate wave functions. 
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Interacting charged particles in the presence of a harmonic confinement appears naturally in a number of physical 
' circumstances. For example, two planar charged particles in the presence of a constant magnetic field, interacting ions 
, in a trap are governed by the combined effect of Coulomb and harmonic potentials. Unlike the Coulomb and harmonic 
cases, the above artificial atom, known as Hooke's atom is not exactly solvable [1-4]. The fact that, a countable infinity 
of eigenstates can be identified exactly for this case, albeit with differing oscillator frequencies, has made Hooke's atom 
a very interesting test ground for checking the effect of correlation arising due to the interaction. The model has been 
used for explicating the working of density functional theory [5]. The recent interest in entangled systems, involving 
continuous variables, makes this model an ideal candidate for a careful study. It should be mentioned that, the 
comparably simpler system of two-interacting oscillators [6] have been recently reanalyzed to see the origin of the 
correlation effects, since in this case both the exact and Hartree-Fock solutions are known [7]. 

After perturbative and numerical investigations of the Hookean atom [1—4], Kais et al., were able to identify one 
state analytically for a specific value of the spring constant [8]. The same system, with an additional linear term 
in the inter-particle potential, was also analyzed by Ghosh et al., to study the single particle density [9]. Later on, 
Taut studied this system carefully and through the analysis of the three term recurrence relation showed that, the 
Q H . system possesses analytical solutions for a particular, infinitely denumerable set of oscillator frequencies [10-12]. He 
further showed that, two conditions need to be satisfied for obtaining the above solutions. The first one relates the 
energy with the oscillator frequency and other quantum numbers of the system. The general expression for the second 
condition, which connects oscillator frequency with the Coulomb coupling and other quantum numbers have not been 
obtained exactly; the solutions for the first few values have been worked out. Recently Holas et al., have studied the 
general structure of the density matrix for this system, concentrating on a single value of the spring constant earlier 
considered by Kais et al [13]. 

In this paper, we provide a general algebraic procedure to find the eigenstates of planar Hookean atom, which can 
be extended to higher dimensions. Our procedure allows one to compute the series expression of the wave function 
explicitly in a more economic manner as compared to the routinely used series solution approach. It is based on a 
novel method for obtaining the solutions of linear differential equations [14, 15]. Apart from obtaining the analytical 
expressions, in the present approach one can also carry out a perturbative expansion for the wave function, when 
closed form solutions are not available. The analytical expressions are then used to compute the exact single particle 
densities for arbitrary values of the quantum numbers. The densities and the information entropy, associated with 
the pair correlation function are systematically studied to analyze the effect of the interaction on correlation in this 
artificial atom. 

In the following section, we outline the above mentioned method for solving linear differential equations and employ 
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the same to systematically study the eigenstates of the Hooke's atom. The single particle densities for the planar 
Hooke's atom for arbitrary values of azimuthal quantum number are obtained analytically. The single particle densities 
and information entropy, associated with the relative motion, are then analyzed. The effect of interaction on correlation 
is shown to be revealed by information entropy. We then point out, how the connection of Hooke's atom with the 
well-studied quasi-exactly solvable (QES) system can be used to develop perturbative expansion involving Coulomb 
coupling parameter Z. We conclude in section III after pointing out directions for future studies. 



II. HOOKE'S ATOM 



The Hooke's atom, dealing with two interacting charged particles in harmonic potential, is governed by the Hamil- 
tonian (in units, K = m = e = 1): 
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This Hamiltonian decouples in the center of mass R = |(ri + r 2 ) and the relative coordinate, r = r 2 — ri, which 
jive rise to the new momentum operators, p = — iV r = \ (p 2 — Pi) and P = — iVr = p 2 + pi- 
Eq.(l) in these coordinates reads 
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where A R = 2A(R) = A(ri) + A(r 2 ), A r = |A(r) = |(A(r 2 ) - A(ri)) and u R = 2w , cj r 
function factorizes: 
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It is clear from Eq.(3) that, the eigenstates of Hr are identical with the two dimensional oscillator states, where 
Ar is a linear function of R, for a constant magnetic field B. 

The Schrodinger equation H r ip(r) = sip(r) for the internal motion, can be cast in the form, 
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Here ujl is the Larmor frequency. The angular and radial part of the wave function y(r), are decoupled through the 
ansatz, 
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where u(r) satisfies the radial equation, 
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Here Lu r = \\Ju\ + lu 2 . 

To solve the above equation, one introduces the dimensionlcss variable p = \fuT r r and 
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and substitute 
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This yields the following equation: 
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The above radial equation can be solved for a series solution, following a recently developed method for solving 
linear differential equations [14, 15]. A single variable differential equation, after suitable manipulations, can be cast 
in the form, 
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here F(D) is a function of the Euler operator D = x-^, possibly containing a constant. P(x,d/dx) contains other 
operators of the differential equation concerned. The series solution for y(x) can be written in the form: 



y{x) = c x \Y i (-\y 



F(D) 



P(x, d/dx) 



(12) 



provided F{D)x x = 0. 

The proof of this can be be checked by direct substitution. For the Hooke's atom case, we have F(D) = D(D + 2m), 
A = 0, -2m and P = (e' - 2\m\ - 2) p 2 - 2p 3 -f - after multiplying p 2 to the Eq.(10). Eq.(12) yields, 
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which can be expanded in the form, 

Z 



tip) = 1 + 



E 



ov(2ra+l) I u r {2m+ l)2(2m + 2) 2(2™ + 2) 
1 f Z 3 ZE 



VS r 1 w r (2m + l)2(2m + 2)3(2m + 3) (2m + l)3(2m + 3) 
ZE 2Z ) , 



2(2m + 2)3(2m + 3) (2m + l)3(2m + 3) 



here E = e' - 2(\m\ + 1) . 



(14) 



In order that, the above series terminates at p n 1 , the coefficients of p n and p n+1 terms must be zero. There are two 
operators |(e' — 2|m| — 2) p 2 — 2/9 3 ^| and -j£=, which increase the degree of the polynomial solution. Once a n , the 

coefficient of p n is zero, p n+1 term will not get any contribution from -^£= and the only contribution to a n+ i will 

arise from the action of |(e' — 2\m\ — 2) p 2 — 2/? 3 ^| on a„_ip™ _1 . Hence, a„+i = yields, e' = 2(n + \m\). The 

coefficients a n can be straightforwardly determined in our approach; although it becomes tedious for higher values of 
n. Once the above condition is implemented, one needs to check if a n — yields exact solutions. 
For n = 2 and arbitrary m and Z, we have the frequencies 
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and 



u(r) = 1 + 



Zr 
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In case of n = 4, we have two roots for the oscillator frequency: 
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For the lower frequency the polynomial part of the wave function is given by, 
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with energy e = 2(n + |m|)a; r . The other eigenvalues and eigenstates can be found analogously. We note that, for 
the general case the coefficients of p n+1 yields a relation between the energy, frequency and the azimuthal quantum 
number. It is the second relationship i.e., a n = which brings the relationship between Z, Co r and m. It is worth 
mentioning that, using the above series expansion, it is possible to approximately determine the other analytically 
inaccessible states of a given Hamiltonian. Although, we intend to deal with this in greater detail later, below we 
briefly outline the underlying procedure. To be specific, consider the repulsive n — 2 case, for which the ground 
state is analytically available for a particular frequency. To obtain approximately, the first excited state of the same 
Hamiltonian, one starts with the above series expansion, retaining sufficient number of terms to ensure that the wave 
function has one real zero in the half line. In this wave function the energy can be treated as a variational parameter, 
to accurately determine the eigenvalue and eigen state, as has been shown for the QES sextic oscillator previously 
[16]. This process will be made more transparent while demonstrating the connection of the Hooke's atom with the 
above mentioned sextic oscillator problem below. 

We now proceed to compute the pair correlation function and single particle density defined respectively as, 
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The expression for the charge density for n = 2, m = and Z = 1 can be computed: 
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We give below the three dimensional plots of information entropy densities associated with the pair correlation 
function G(r) = |</?(r)| 2 , as the same incorporates the interaction, defined as Sg(y) = — G(r) ln[G(r)], for various 
values of n and Coulomb couplings. The effect of interaction on the correlation is clearly seen in these plots. In 
Fig.l we have plotted Sg{?) for n = 2, Z = —1 and various values of azimuthal quantum number m = 0, 1 and 2. 
We observe a dip at the origin; this occurs because G(r) > 1 at r = which makes the quantity 5c(r) < 0. This 
is a special case, as for higher n values and m = 0, we observe a peak at the origin, since G(r) < 1. For higher 
values of m correlation decreases and wave function of the system delocalizes. Similar behavior is seen in Fig. 2. Fig. 
3 explicates the effect of Coulomb interaction on the entropy density; Fig. 3 (a), (b) and (c) are plotted for various 
coupling strengths keeping n and to fixed. It is clear that, as the coupling strength increases, the system becomes 
more localized. In Fig. 4, we have plotted the entropy density for repulsive Hooke's atom with n = 4. We observe, 
for the m = there is a small dip at the center, which becomes wider for higher values of to, which agrees with the 
fact that for higher values of to there is less correlation present in the system. 

In Fig. 5, we have numerically calculated the position space entropy for various values of to, keeping n — 3, for 
attractive and repulsive cases. It is clear from the plot that, as the oscillator frequency increases (which corresponds 
to a smaller value of to), entropy decreases. Also an attractive system has more entropy as compared to its repulsive 
counterpart. 

The Hooke's atom can be connected to the quasi-exactly solvable models by a simple coordinate transformation. 
We can make use of this connection to develop a suitable perturbation theory involving the coupling parameter Z. 
The QES systems are intermediate to exactly solvable and non-solvable systems. Unlike the exactly solvable models, 
in such systems only a finite number of states can be determined analytically [17, 18]. The typical example of a 
QES system, which we are going to map with Hooke's atom is a sextic oscillator with a centrifugal barrier, whose 
Hamiltonian is given by, 
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We make the following similarity transformation H = ip l (x)Hi[>o(x), with ^0(2;) = x m+1 cxp(— x^y/j/A) and seek 
the polynomial solution for the reduced Hamiltonian, 
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The series solution of the above equation, in the present approach, can be written for F(D) = D(D + 2m + 1): 
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FIG. 1: Plots of information entropy density for attractive Z = —1 planar Hooke's atom with n — 2 and (a) m = 0, (b) m = 1 
and (c) m = 2. 

where, ,4 = (§ + ^ + (to + 1)V7) • 

It is worth noting that, the operators in the above equation, y/^x 3 -^ and + + (to + 1)^/7^ a; 2 increase the 

degree of the polynomial by two. In order to preserve the degree of the polynomial in an eigenvalue equation of H, 
the following condition on the parameters should be imposed, 
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A simple transformation of variable: x 2 = r in Eq.(24) allows us to map this QES problem to Hooke's atom; 
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With ^/;(r) = r 1,/4 </?(r) and (to 2 — 1/4) = m ( r " +1 ) — JL the above equation, comes to the form of Hooke's atom, 
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It is interesting to notice that, energy eigenvalue, quadratic and sextic couplings of the QES Hamiltonian become 
the Coulomb coupling, energy and quadratic coupling of the corresponding Hookean system respectively, but the 
centrifugal barrier term remains same in both the cases. With this parameterization, it can be easily seen that the 
QES condition given by Eq.(28) is same as the reduced energy e' = 2(n + |m|). As has been mentioned earlier, for 
finding the analytically inaccessible QES states, one can start with the above series expansion and terminate it at a 
desired degree, depending on the required number of the nodes of the concerned wave function. Treating energy as a 
variational parameter, one minimizes {^\(H — E)\ty), to obtain an accurate expression for the eigen state and eigen 
value [16]. The perturbation in energy in QES regime implies a perturbation in Z in its Hookean counterpart. 



III. CONCLUSIONS 



In conclusion, we have illustrated the utility of a recently developed method for solving linear differential equations 
to Hooke's atom The information entropy and their densities are analyzed systematically for studying the effect of 
interaction on correlation. The procedure for developing perturbative expansions based on the present approach is 
indicated. A connection of this dynamical system with well studied QES systems can also be used for developing a 
suitable perturbation theory involving the charge parameter Z. The usefulness of single particle density to density 
functional theory is well-known, as also its ability for studying entanglement of this correlated system. We would like 
to get back to these questions in near future. Recently, Ralko and Truong [19] has connected this system to anyons 
and Heun's equation [20]. This is an interesting direction, needing further investigations. 
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FIG. 2: Plots of information entropy density for attractive Z = — 1 planar Hooke's atom with n = 3 and (a) m = 0, (b) m = 1 
and (c) m = 2. 




FIG. 3: Plots of information entropy density for attractive planar Hooke's atom with n = 2 and m = depicting the effect of 
Coulomb coupling strength for (a) Z — —2, (b) Z = — 3 and (c) Z = —7. 




FIG. 4: Plots of information entropy density for repulsive (Z = 1) planar Hooke's atom with n = 4 (a) m = 0, (b) m = 1 and 
(c) m = 5. 
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FIG. 5: Plot depicting the variation of position space entropy with respect to various allowed oscillator frequencies for n = 3, 
the green line and red lines are for attractive and repulsive cases respectively. 
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